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We report on the inclusive branching fractions of and of B'^ mesons decaying to D^X, D^X, 
D'^ X, X, Dfx, X, Af X, X, based on a sample of 88.9 million BB events recorded with 
the BABAR detector at the T(45) resonance. Events are selected by completely reconstructing one 
B and searching for a reconstructed charmed particle in the rest of the event. We measure the 
number of charmed and of anti-charmed particles per B decay and derive the total charm yield per 
B- decay, n" = 1.313 ± 0.037 ± 0.062l^:2|^, and per B° decay, n° = 1.276 ± 0.062 ± Q.Qh&toMQ 
where the first uncertainty is statistical, the second is systematic, and the third reflects the charm 
branching-fraction uncertainties. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



The dominant process for the decay of a quark is 
h cW*~ [1], resulting in a (flavor) correlated c quark 
and a virtual W. In the decay of the W, the produc- 
tion of a or a cs pair are both Cabibbo-allowed and 
should be equal, the latter being only suppressed by a 
phase-space factor. The first process dominates hadronic 
h decays, while the second can be easily distinguished as 
it will produce a (flavor) anti-correlated c quark. Exper- 
imentally, correlated and anti-correlated charm produc- 
tion can be investigated through the measurement of the 
inclusive i?-decay rates to flavor-tagged charmed mesons 
or baryons. Current measurements [2-4] of these rates 
have statistically limited precision and do not distinguish 
among the different B parent states. 

Most of the charged and neutral D mesons produced 

in B decays come from correlated production B DX. 
However, a significant number oi B ^ DX decays are 
expected through h — > ccs transitions, such as B ^ 
£)(*) X(*\mT) . Although the branching fractions of 
the 3-body decays B D^*^D^*''K have been mea- 
sured [5, 6], it is not clear whether they saturate B — > 
DX transitions. It is therefore important to improve the 
precision on the branching fraction B{B DX). 

By contrast, the anti-correlated Dj production 
B — > D~D{mr) is expected to dominate B decays to 
mesons, since correlated production needs an extra ss 
pair created from the vacuum to give B — > DfK~{mr). 
There is no prior published measurement of B{B — > 
D+X). 

All strangeless charmed baryons decay to Ac- Corre- 
lated Ac are produced in decays like B~ — > yl+p7r~(7r), 
while anti-correlated A~ should originate from B~ — > 
ScA~{Tr). Another possibility is B~ AfA~K~, the 
baryonic analogue of the DDK decay. The rates for Ec 
production in B decays [7] are unknown, because there 
is no absolute measurement of Sc decay branching frac- 
tions. 

This analysis uses T{4S) BB events in which either 
a -B+ or a B° meson (hereafter denoted B^^co) decays 
into a hadronic final state and is fully reconstructed. We 
then reconstruct D, Ds and Ac from the recoiling B^ 
(i?°) meson and compare the flavor of the charm hadron 
with that of the Bieco, thus allowing separate measure- 
ments of tlie B- (_B") ^ D° X, D+ X, p+ X, A+ X 
and B- (5°) ^ D° X, X, X, A- X branch- 



ing fractions. We extract B{B~ AfAcK^) from the 
missing-mass spectra of the AfK^ or Ac systems re- 
coiling against the B^eco- We can then evaluate indirectly 
B{B- ScX) = B{B- -> A-X) - B{B- A+A-R-) 
and compute the average number of charm (anti-charm) 
particles per B~ decay, N~ (N.^): 

N- = Y^BiB-^XcX), (1) 
= Y.^{B-^XcX), (2) 

where the sum is performed over Xc = D^, D^ , _D+, 
yl+, Sc, (cc) or Xc = , D°, Dj , Ac, (cc), and (cc) 
refers to all charmonium states collectively. We neglect 
Ec production, as it requires both a cs and an ss pair in 
the decay to give Sci^c- We can sum and to ob- 
tain the average number of charm plus anti-cliarm quarks 
per B~ decay, n~ = N~ + N.^ (and similarly for B° de- 
cays). In addition to the theoretical interest [8-10], the 
fact that anti-correlated charmed particles are a back- 
ground for many studies also motivates a more precise 
measurement of their production rates in B decays. 

The measurements presented here are based on a sam- 
ple of 88.9 million BB pairs (81.9 fb"^) recorded at 
the T{AS) resonance with the BABAR detector at the 
PEP-II asymmetric-energy B-meson factory at SLAC. 
The BABAR detector is described in detail elsewhere [11]. 
Charged-particle trajectories are measured by a 5-layer 
double-sided silicon vertex tracker and a 40-laycr drift 
chamber, both operating in a 1.5-T solenoidal magnetic 
field. Charged-particle identification is provided by the 
average energy loss (dE/dx) in the tracking devices and 
by an internally reflecting ring-imaging Cherenkov de- 
tector. Photons are detected by a CsI(Tl) electromag- 
netic calorimeter. Wc use Monte Carlo simulations of 
the BABAR detector based on GEANT4 [12] to optimize 
selection criteria and determine selection efficiencies. 

We reconstruct B"*" and B" decays {B^cco) in the modes 
B+ ^ D(*)07r+, 5(*)V, -D(*)Oa+ and B^ ^ D^*)-tt+, 
D<^*'^-p+, D^*^-at. -D° candidates are reconstructed in 
the K+TT-, K+n-n", K+t:-tt+tt- and K'^tt+ tt' {K'^ -> 
7r+7r~) decay channels, while D~ are reconstructed in 
the Jsr+7r~7r~ and K'^'k~ modes. D* candidates are re- 
constructed in the D*' D^tt- and _D*" _D°7r", 
£)°7 decay modes. The first kinematic variable used to 
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identify fully reconstructed B decays is the beam-energy 
substituted mass, rriEs = a/ (s/2 + Pj • PB)^/Ef — p^, 
where pb is the -Breco momentum and (Ei,pi) is the 
four-momentum of the initial e+e~ system, both mea- 
sured in the laboratory frame. The invariant mass of 
the initial e+e^ system is ^/s. The second variable 
is AE = Eg — where Eg is the Breco candi- 

date energy in the center-of-mass frame. We require 
\AE\ < naAE with n = 2 or 3, depending on the de- 
cay mode, and using the measured resolution (tae for 
each decay mode. 

In the TTiES spectra (Fig. 1), we define a signal re- 
gion with 5.274 < tties < 5.290 GeV/c^ and a back- 
ground control region with 5.220 < mEs < 5.260 GeV/c^. 
For each of the B-decay modes, the combinatorial back- 
ground in the signal region is derived from a fit to 
the Toes distribution that uses an empirical phase-space 
threshold function [13] for the background, together with 
a signal function [14] peaked at the B meson mass. The 
numbers of reconstructed B^ and B^ candidates, Ng+ — 
85840 ± 1910 (syst.) and Ngo = 48322 ± 590 (syst.), are 
then obtained by subtracting this background from the 
total number of events found in the signal region. These 
measured B meson yields provide the normalization of all 
branching fraction measurements reported below. The 
systematic uncertainties quoted above are computed by 
varying the boundaries of the signal and background re- 
gions, and by comparing the shapes of the threshold func- 
tion [13] in the data and in the simulation. 

The contamination of B^ events in the i?+ signal in- 
duces a background which peaks near the B mass. From 
the Monte Carlo simulation, the fraction of 5° events 
in the reconstructed B~^ signal sample is found to be 
Co — 0.034, and the fraction of B'^ events in the recon- 
structed B° signal sample to be c+ = 0.019. A 100 % 
systematic uncertainty is conservatively assigned to these 
numbers but they will have a small effect on the final re- 
sults. 

We now turn to the analysis of inclusive D, Dg and Ac 
production in the decays of the S's that recoil against the 
reconstructed B. Charmed particles Xc (correlated pro- 
duction) are distinguished from anti-charmed particles 
Xc (anti-correlated production). They are reconstructed 
from charged tracks that do not belong to the Sreco- The 
decay modes considered are listed in Table I. 

For charged B decays. Fig. 2 shows the D, Dg, and 
Ac mass spectra of correlated and anti-correlated candi- 
dates recoiling against B's reconstructed in the toes sig- 
nal region, for some selected decay modes. These spec- 
tra are fitted with the sum of a Gaussian signal and a 
linear background (including a satellite peak for some 
channels [15]). The shaded areas correspond to well re- 
constructed D, Ds or Ac from the combinatorial B-ceco 
background. They are obtained from data in the toes 
background control region, normalized to the number 
of combinatorial background events expected under the 
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FIG. 1: mBS spectra of reconstructed (a) B'^ and (b) 
candidates. The full vertical line shows the upper limit of 
the background control region (hatched), the dotted vertical 
line the lower limit of the B signal region. The crossed area 
shows the background under the B signal. The solid curve is 
the sum of the fitted signal and background, the dashed curve 
is the background component only. 



Breco peak. The background-subtracted reconstructed 
signal yields are listed in Table I. The reconstruction 
efficiencies for each charmed (anti-charmed) final state 
Xc — + / {Xc — *■ /) are computed from the simula- 
tion as a function of the charmed-particle momentum 
in the B^ center-of-mass frame, and are applied event- 
by-event to obtain the efficiency-corrected charm signal 
yields N{Xc /) {N{Xc /)). The final branch- 
ing fractions are computed from these yields, the num- 
ber of Bj-cco, and the intermediate branching fractions 
B{Xc /) taken from [16]. They are given by 



BiB- ^ XcX) = 



NjXc ^ /) 
Nb+ X B{Xc ^ ,/) 



-coBo. (3) 



Here the raw branching fraction for B^ XcX is mod- 
ified by a small corrective term, cqBq, that accounts for 
the i?^ contamination in the reconstructed B~^ sample. 
The factor Bq depends on the measured B^ XcX and 
B'^ XcX branching fractions, and on the B'^ — B^ mix- 
ing parameter Xd [16]. It ranges from less than 3% for 
Ac to as much as 50% for correlated D° and Z?"*". Dou- 
bly Cabibbo-suppressed decays are also taken into 
account. The branching fractions and their errors are 
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TABLE I: Charmed-particle signal yields and B branching fractions per decay mode. The first uncertainty is statistical, the 
second is systematic (but does not include the charm branching fraction uncertainties). 





















Xc decay mode 






ts 








dO 

D 


V V 

— > AcA 




yield 


B{%) 


yield 


B{%) 


yield 




yield 






1273±42 


79.2±2.6±3.9 


160±16 


9.3±1.0±0.5 


397±24 50.3±3.4±2.4 


139±14 


7.3±2.2±0.5 




998±65 


80.6±5.3±7.5 


173±30 


13.4±2.4±1.3 


332±36 


56.2±6.8±5.4 


83±23 


1.8±4.4±0.5 


D+^K-n+n+ 


262±29 


9.8±1.2±1.2 


98±20 


3.8±0.9±0.4 452±31 


39.7±3.0±2.8 


125±18 


2.3±1.8±0.3 




11±5 


2.2±1.1±0.3 


82±11 


16.5±2.3±1.7 


24±6 


8.3±2.8±0.8 


28±6 


9.9±2.9±1.0 




0±3 


0.0±1.1±0.2 


55±11 


18.0±3.5±1.7 


3±4 


0.0±2.8±0.1 


14±5 


9.9±4.1±1.2 




0±3 


0.0±0.9±0.2 


31±9 


9.2±2.7±0.8 


12±5 


5.0±3.4±0.4 


23±6 


13.3±4.3±1.0 




41±9 


3.5±0.8±0.3 


33±9 


2.9±0.8±0.3 


28±8 


4.9±1.7±0.4 


16±6 


2.0±1.2±0.2 



given in Table I. The statistical and systematic uncer- 
tainties are computed separately for each channel. For 
example, the 3.9% absolute systematic uncertainty on 
B{B~ D'^{K^'!t^)X) reflects the quadratic sum of 
1.8% attributed to Nq-i-, 1.3% to the error on the rate 
of true D's in the B combinatorial background, 0.8% to 
the Monte Carlo statistics, 1.2% to the track- finding effi- 
ciency, 2.5% to the particle identification, 1.2% to cq, and 
0.1% to Bq. We combine the results from the different 
£)" and Ds decay modes to extract the final branching 
fractions listed in Table II. 



TABLE II: Combined B~ branching fractions. The first un- 
certainty is statistical, the second is systematic, and the third 
reflects charm branching- fraction uncertainties [16]. 





correlated 


anti-correlated 




B{B- ^ X,X){%) 


B{B- ^XcX)[%) 


D° 


79.3 ± 2.5 ±4.0ti:9 


9.8 ± 0.9 ± 0.5lo;-3 


D+ 


9.8 ± 1.2 ± 1.2to;7 


3.8±0.9±0.4i°;^ 


Dt 


0.5 ± 0.6 ± Q.2toX 
< 2.2 at 90% CL 


14.3 ± 1.6 ± 1.5l3:^ 


At 


3.5±0.8±0.3lJ:| 


2.9 ± 0.8 ± 0.3to.6 



To extract Nc from these numbers, we need to evalu- 
ate the contribution of . Combining the 
four-momenta of the recoiling B~ , of a K~ and of the 
reconstructed Af or A~ candidate, we compute the miss- 
ing mass: the absence of signal at the Ac mass excludes a 
significant contribution of this process. We therefore take 
B{B- ScX) = B{B- A'X) in the computation of 
Nc- Using Eqs. 1 and 2 and taking B{B~ — > {cc)X) = 
(2.3 ±0.3)% [17] [18], one obtains: 

A^- = 0.983 ± 0.030 ±0.046to 023, 
AT- = 0.330 ± 0.022 ± 0.020l°:^^J, 
n- = 1.313 ±0.037 ±0.0621°;°^^. 

The reconstruction of D, Ds and Ac from B° decays is 
performed in the same way as that in the B~ analysis. 



The corresponding yields are listed in Table I. We then 
compute for each decay channel Xc f the efficiency- 
corrected signal yields N{Xc /) {N{Xc /)) and 
define the raw branching fractions Be and Be as 

Be = N{Xc ^ f)/[NBO X B{Xe ^ /)] (4) 
Be = N{Xe ^ f)/[Nso X BiXe ^ /)]. (5) 

After correcting these numbers for B'^B^ mixing, we 
obtain the final branching fraction for B*^ — > XeX: 

B{B'^XeX) = Bo-xdBc + Bc)-c+B+ ^ 

1 - 

where Xd = 0.181 ± 0.004 is the B" - _B° mixing pa- 
rameter [16]. The correcting factor B+ accounts for 
5+ contamination in the B*^ sample and depends on 
B{B~ XcX) and BiB+ XcX). The results are 
given in Table I. Combining the different 1)° or 
modes, we obtain the final branching fractions listed in 
Table III. 



TABLE III: Combined B° branching fractions. The first un- 
certainty is statistical, the second is systematic, and the third 
refiects charm branching-fraction uncertainties [16]. 





correlated 


anti-corrclatcd 


Xc 


BiB" ^ XcX){%) 


B(B" ^ XcX){%) 


D° 


51.1 ±3.1 ±2.5ljj 


6.3 ± 1.9 ± 0.51^ 2 


D+ 


39.7 ± 3.0 ± 2.8tli 


2.3 ± 1.8 ± 0.31^ 2 

< 5.1 at 90% CL 


Dt 


3.9± 1.7±0.4tJ;| 
< 8.7 at 90% CL 


10.9 ± 2.1 ±o.8l^:| 


At 


4.9 ±1.7 ±0.411* 


2.0 ± 1.2 ± 0.21°;;^ 
< 3.8 at 90% CL 



To compute A^c, we neglect B^ — > AtA^ produc- 
tion and assume that B{B^ EcX) = B{B^ A'X). 
Substituting B^ for B~ in Eqs. 1 and 2 and taking 
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FIG. 2: Correlated (left) and anti-correlated (right) charmed 
particle mass spectra in the recoil of B'^ events, for (a),(b) 
D° K-n+- (c),(d) D+ K--K+-K+; (e),(f) D+ (j}Tv+; 
and (g),(h) A'^ — > pK~it^ . The solid curve is the sum of 
a Gaussian signal and of a linear background plus mode- 
dependent satellite contributions [15]. The shaded areas show 
the contribution of well reconstructed D, Ds or Ac in the 
combinatorial background. 



in Table IV. We obtain an upper limit on the correlated 
L>+ fraction in B' decays : B{B- D+X)/B{B- 
Dfx) < 0.126 at 90% CL. 

TABLE IV: Fraction w of anti-correlated charm. 



Mode 


B decays 


B''^ decays 


D'^X 


0.110 ±0.010 ±0.003 


0.110 ±0.031 ±0.008 


D-X 


0.278 ± 0.052 ± 0.009 


0.055 ± 0.040 ± 0.006 


D-X 


0.966 ± 0.039 ± 0.012 


0.733 ± 0.092 ± 0.010 


A-X 


0.452 ± 0.090 ± 0.003 


0.286 ±0.142 ±0.007 



In conclusion, we have measured for the first time the 
branching fractions for inclusive decays of B mesons to 
flavor-tagged Z?, Dg and separately for B~ and B^ . 
We observe significant production of anti-correlated 
and _D+ mesons in B decays (Table IV) , with the branch- 
ing fractions detailed in Tables II and III. The correlated 
Ds production in B~ decays is measured to be small. 

As expected, the sum of all correlated charm branch- 
ing fractions, A^c, is compatible with 1, for charged as 
well as for neutral i3's. The numbers of charmed parti- 
cles per B- decay (n" = 1. 313 ± 0.037 ±0.062^!^;^^^) and 



1.276 ±0.062 ±0.058 



+0.066\ 
^0.046-' 



per B^ decay (r 
sistent with previous measurements [2, 17, 19] and with 
theoretical expectations [8-10]. 
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B{B° (cc)X) = (2.3 ± 0.3)% [17] [18], we obtain: 

7V° = 1.039 ± 0.051 ±0.049l^;^i^?, 
= 0.237 ± 0.036 ±0.012i;] J]^^, 



n° = 1.276 ±0.062 ±0.058 



+0.066 
-0.046- 



We also compute the fraction of anti-correlated 
charm production in B decays, 'w{Xc) = B{B 
XcX)/{B{B XcX) ± B{B XcX)). Here, many sys- 
tematic uncertainties cancel (tracking, K identification, 
D branching fractions, B counting). The results are given 
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